Transport in Ultraclean YBa2Cu307: neither Unitary nor Born Impurity Scattering 
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The thermal conductivity of ultraclean YBa2Cu307 was measured at very low temperature in 
magnetic fields up to 13 T. The temperature and field dependence of the electronic heat conductivity 
show that two widespread assumptions of transport theory applied to unconventional superconduc- 
tors fail for clean cuprates: impurity scattering cannot be treated in the usual unitary limit (nor 
indeed in the Born limit), and scattering of quasiparticles off vortices cannot be neglected. Our 
study also sheds light on the long-standing puzzle of a sudden onset of a "plateau" in the thermal 
conductivity of Bi-2212 versus field. 

PACS numbers: 72.15.Eb, 74.72. Bk, 74.25.Fy 



Since 1985, the theory of transport in unconventional 
superconductors has been dominated by the ubiquitous 
assumption that impurity scatteringmust be treated in 
the unitary limit (7r/2 phase shift) The aim of this 
Letter is to test this assumption in the simplest and most 
reliable context available, namely that of the cuprate 
material YBa2Cu30a; (YBCO), for which the supercon- 
ducting gap structure is firmly established and straight- 
forward, with d.j.2_y2 symmetry over an approximately 
cylindrical Fermi surface. We do this by measuring heat 
transport in ultraclean single crystals that have a scatter- 
ing rate one order of magnitude lower than in previous 
studies. This allows us to reliably resolve the intrinsic 
temperature dependence of the electronic thermal con- 
ductivity. 

The pivotal result of transport theory is the universal 
limit of conductivity, whereby the ability of a supercon- 
ductor to carry heat as T — > is independent of either 
the concentration of scattering centers or the strength of 
the scattering potential {i.e. phase shift). The univer- 
sal limit was first observed in YBCO 2], and then con- 
firmed in Bi2Sr2CaCu208 (Bi-2212} Q, as well as in the 
p-wave superconductor Sr2Ru04 It is only by going 
beyond the universal limit that one can test the assump- 
tion of imitary scattering. The increase in conductivity 
with temperature T or applied magnetic field H is ex- 
ected to depend strongly on scattering phase shift 6 
The main finding of the present study is that the 
electronic thermal conductivity increases with T much 
faster than could ever be expected from unitary scatter- 
ing. The high degree of order in the present crystal also 
reveals unambiguously the dominance of vortex scatter- 
ing over impurity scattering even at modest fields. This 
will shed light on the long-standing puzzle of a sudden 
onset of a "plateau" observed previously in the thermal 
conductivity of Bi-2212 versus field 7, 8]. 

The thermal conductivity k was measured using a sin- 
gle heater-two thermometer method. The heat current 



was supplied along the a-axis of the sample, and the 
magnetic field applied parallel to the c-axis. The mea- 
surements were made in a dilution refrigerator by vary- 
ing the temperature from 0.04 K to > 0.7 K at fixed 
magnetic field. The samples were field-cooled by cycling 
to T > 100 K before changing the field. The error in 
the absolute value of the conductivity is estimated to be 
approximately 10%. The relative error between temper- 
ature sweeps at different fields is of order 1%. 

The single crystal platelet of YBCO used in this study 
has dimensions 1.0 x 0.5 mm^ and 25 /ini thick. It was 
grown in a BaZrOa (BZO) crucible 9], which results 
in crystals with extremely high chemical purity (99.99 
- 99.995%) and a high degree of crystalline perfection as 
compared with crystals grown in Y203-stabilised Zr02 
(YSZ) crucibles. The sample was detwinned at 250°C 
under uniaxial stress, followed by a 50 day annealing pro- 
cess at 350°C, resulting in oxygen chains with less than 
0.7% vacancies. This level of oxygen doping, x — 6.99, 
is slightly above that for maximal Tc (93 K), resulting in 
a marginally lower Tc = 89 K. Four electrical contacts 
were applied using silver epoxy fired at 350°C for 1 hour, 
giving typical contact resistances of 200 m^l at 4.2 K. 

The high purity of our sample can be established from 
its thermal conductivity at high temperature. A low de- 
fect level leads to a larger peak in k below Tc, as noted 
by Zhang et al 10] . In the inset to Fig. we show the 
thermal conductivity for this sample as compared to a 
detwinned YBCO sample with x = 6.95 grown in an 
YSZ crucible. The approximate doubling in peak height 
can be attributed to an order of magnitude decrease in 
the intrinsic disorder level (see also (l"3|). 

The main panel of Fig. ^ shows the temperature de- 
pendence of the thermal conductivity of the high-purity 
(BZO) sample at zero field and in applied fields of 0.8 and 
13 T. The thermal conductivity is the sum of electronic 
and phononic contributions: k = k,, + Kph- The resid- 
ual linear term, kq/T, obtained by extrapolating k/T to 
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FIG. 1: Thermal conductivity divided by temperature versus 
for a high-purity YBCO sample for magnetic fields applied 
perpendicular to the ab plane. The dashed line is a guide to 
the eye. Inset: High-temperature data for YBCO samples 
grown in different crucibles. 



T — > 0, is entirely electronic. The question is how to ex- 
tract any T-dependence of Ke/T. We use a magnetic field 
to do this. First, note that the zero-field curve in Fig. 
shows a more rapid increase with temperature than the 
two in-field curves which are approximately parallel. On 
the assumption that the phonon transport at very low 
temperatures is limited by scattering from the bound- 
aries of the sample (see and is unaffected by scat- 
tering from vortices, we can only attribute this difference 
to electrons. Furthermore, since all subsequent in-field 
curves lie parallel we assume that this additional elec- 
tronic conduction is completely suppressed when a mag- 
netic field is applied (see inset of Fig.[2Jl, and the remain- 
ing temperature dependence of k/T is due entirely to 
phonons. In other words Kph is the T-dependent part of 
the 13 T data: Kph/T = k(13 T)/T- ko(13 T)/r, where 
Ko(13 T)/r = 0.31 mW/K^cm. The electronic conduc- 
tivity, Ke/T, is then given by subtracting this from the 
total conductivity: Ke{H,T)/T = k{H,T)/T - Kph/T. 
This is shown in Fig.|21for applied magnetic fields from 
0-13 T. 

The zero-field electronic conductivity shows a rapid 
growth with temperature, increasing by a factor of five 
within 0.5 K. The inset plots the same data on a T^ tem- 
perature scale showing that the growth is cubic in tem- 
perature. As soon as a magnetic field is applied this tem- 
perature dependence is completely suppressed (see inset 
of Fig. |2Jl, an effect which can only be attributed to the 
scattering of quasiparticles by vortices. The possibility 
that the change in temperature dependence between zero 
and applied field is due to scattering of phonons by vor- 
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FIG. 2: Electronic thermal conductivity, plotted as k^/T 
vs r, for magnetic fields from zero to 13 T. The electronic 
contribution is extracted as discussed in the main text. Inset: 
The low field curves plotted against T^. 



tices is ruled out by the lack of field dependence above 
0.8 T. By increasing the field to 13 T, an order of magni- 
tude more vortices have been introduced to the system, 
yet the total conductivity remains essentially unchanged. 

Zero magnetic field (H—0). Using a self-consistent 
quasiclassical theory, formulated at low temperatures 
where heat transport is limited by electron scattering 
from random defects [isl , the thermal conductivity 
due to quasiparticles at the nodes of a d-wave supercon- 
ductor is given by 



1l (r) ^ 1^ 
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where kqo/T is the universal conductivity limit, 7 is the 
impurity bandwidth and the coefficient a is strongly de- 
pendent on the scattering phase shift. This expression 
is valid in the dirty limit where ksT < 7. The tem- 
perature dependence of our extracted electronic conduc- 
tivity is well described by this form. Fitting this ex- 
pression to the zero field data (see inset of Fig. ^ gives 
Ke/T = 0.16(1 + 19.2T^). Note that this corresponds to 
a huge 20-fold increase in Ke/T by 1 K. 

H = 0: Zero temperature. The residual linear term, 
kq/T = 0.16 mW/K^cm, is in excellent agreement 
with the value published for optimally doped YBCO 
(0.14 mW/K^cm ^). This is direct confirmation of the 
universal nature of low temperature quasiparticle thermal 
conductivity. In this case we have measured a sample in 
the more difficult regime of increased purity, where we 
observe an order of magnitude lower scattering rate than 
previously, and still recover the universal limit. (A slight 
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increase in kq/T is expected based on the small increase 
in doping from x = 6.95 to x = 6.99 [l^.) 

H = 0: Finite temperature. From the temperature 
dependent part of our fitted data we can estimate the 
impurity bandwidth and scattering rate. For scattering 
in the unitary limit a = 1/2 and the impurity bandwidth 
7 is related to the normal-state scattering rate r„ by 
the relation 7 = 0.63^ AoTn Using the value from 

the fitted zero-field data (inset of Fig.|3), we obtain 7 ^ 
0.25 K and using Aq = 2.14 ksT^, we get r^/T^ ~ 10-^ 
therefore r„ ~ 1 x 10^ s~^. Such a small scattering 
rate is unrealistic. Using vp — 2.5 x 10^ cm/s 12], it 
would imply a normal-state mean free path as long as 
the longest dimension of the sample: I ^ 1 mm! 

In the Born approximation, a — (ttv2To)/2 and 7 = 
4Ao exp(— 7rAo/2r„). Assuming a pure d-wave gap gives 
V2 = 2Ao/hkp and using tq = l/2r„, leads to 7 ~ 3 K 
and r„ ~ O.6A0 ~ 2.5 x 10^^ s^^. Again we estimate a 
scattering rate that is unrealistic, in this case much too 
large. If the scattering rate were truly this magnitude it 
would lead to a substantial suppression of Tc, as noted 
previously [l^ . which is not observed experimentally. 

From the measured temperature dependence of the elec- 
tronic thermal conductivity, we conclude that the usual 
quasiclassical calculation cannot be correct if it treats im- 
purity scattering with a single isotropic phase shift of ei- 
ther (Born) or 7^/2 (unitary). 

In a broader context, similar measurements on uncon- 
ventional superconductors UPta and Sr2Ru04 re- 
veal conductivities far too small to consider Born scat- 
tering. They are therefore analysed in the unitary limit. 
Consistent with the present work, the trend that emerges 
is that a quantitative analysis of transport in the super- 
conducting state, using the unitary limit, leads to mean 
free paths considerably longer than are consistent with 
normal state measurements such as de Haas van Alphen 
studies or resistivity above Hc2. This suggests that the 
standard theoretical approach to transport in unconven- 
tional superconductors is generically inadequate. 

Microwave conductivity measurements |l6lll7| on sam- 
ples of identical quality to that used in this study have 
also been compared to current theories of transport in un- 
conventional superconductors [l^ . These measurements 
refiect some but not all of the characteristics of weak 
scattering (Born limit) and point to either inadequacies 
in the conventional theories, or the need to consider in- 
termediate phase shifts, or both. 

An earlier study on Bi-2212 reported an electronic 
conductivity, Ke/T, that increased as T rather than 
[3 . The phonon contribution is subtracted by comparing 
measurements before and after electron irradiation. A 
linear T-dependence is expected in the clean limit, where 
7 < ksT. It is unclear, however, why Bi-2212 crystals 
with a level of disorder two orders of magnitude larger 
than our BZO-grown YBCO crystals should be in that 
limit. Indeed it is a surprise that any T-dependence of 
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FIG. 3: Normalized residual electronic conductivity as a func- 
tion of magnetic field for a slightly overdoped {x = 6.99), 
very pure, detwinned sample of YBCO (triangles). For 
comparison, the same is plotted for optimally-doped YBCO 
{x — 6.95), where the sample is less pure 19], along with 
the fit to a semiclassical theory which describes the data 
(dashed line). An attempt to fit the same theory to the BZO- 
grown sample is also shown (dotted line). Inset: High tem- 
perature isotherm at 0.5 K for the pure sample compared to 
a magnified (x 40) plot of an isotherm at 6 K for Bi-2212 0. 

Ke/T can be resolved in such crystals. (A 2-fold increase 
is observed by 1 K). 

Finite magnetic field (H > 0). In Fig. |31the extrapo- 
lated linear thermal conductivity at T — > is plotted as 
a funtion of magnetic field. The values are normalised 
by the zero-field value kq/T. Also plotted is the previous 
data for an optimally-doped YBCO {x = 6.95) sample 
grown in a YSZ crucible jl9(. In contrast to the lat- 
ter sample, the much purer one shows a rapid increase 
at fields below H = 0.4 T followed by a sudden change 
to a regime where its behaviour is almost field indepen- 
dent. Concomittantly, the coefficient of the term, and 
consequently the quasiparticle lifetime, collapses (see in- 
set of Fig. 121). We interpret this as an indication that 
quasiparticle- vortex scattering is strongly infiuencing the 
transport in this very clean material. 

Current semi-classical treatments of the effect of mag- 
netic field on the low temperature electronic thermal con- 
ductivity predict smooth sublinear increases to within 
logarithmic corrections 0, . This arises when quasi- 
particle energies are doppler shifted by the superfiuid 
flow around magnetic vortices. Earlier measurements, 
reproduced in Fig. |31 on YSZ-grown samples are 
well described by such theories and have been viewed 
as additional evidence for d-wave symmetry of the su- 
perconducting gap. Clearly the field dependence of the 
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measurements reported here on a BZO-grown sample is 
of completely different character and cannot be described 
by available theories. To emphasize this the dotted line 
in Fig. O is the best fit for the theory of Kiibert and 
Hirschfeld jij to the initial rise of the conductivity. In- 
triguingly the scattering rate deduced from this fit is an 
order of magnitude smaller than for the YSZ sample in 
agreement with the relative change inferred from the rise 
in thermal conductivity below Tc. This suggests that 
the rapid growth at low-fields, due to the low impurity 
scattering rate, is truncated by an additional scattering 
mechanism excluded from current theories. 

A common feature of these theories is that at low tem- 
peratures the eS'ect of scattering from vortices is expected 
to be negligible in comparison to impurity scattering. 
However given the demonstrably low impurity scatter- 
ing rate in this BZO-grown sample we argue that such a 
simplification may be incorrect. Assuming these scatter- 
ing mechanisms to be independent, one could add them 
in a Matthiessen-like manner. At zero field the conduc- 
tivity will be in the impurity dominated regime, whilst 
in the very high field limit, vortices will dominate the 
scattering. The crossover between these two limits will 
depend on the relative amount of impurity scattering and 
the cross-section for vortex scattering. Such a model may 
also account for the similar magnitude of the conductiv- 
ity at high fields for the two samples shown in Fig. 
despite their obviously different impurity concentrations. 

A similar strategy was adopted by Franz liH to in- 
clude scattering from a random vortex lattice. In this 
case the theory was developed to explain a plateau-ing 
of the field dependence of the thermal conductivity at 
higher temperatures and has yet to be extended to low 
enough temperatures for comparison with the data here. 
Nevertheless, the phenomenology may be correct; a 
increase coming from a doppler shift is exactly compen- 
sated by scattering from vortices, where the scattering 
length is given by the inter-vortex separation ^ \fH. 
What remains unusual is the sharpness with which this 
compensation turns on. 

With this is mind, the inset to Fig. O shows the field 
dependence of an isotherm at 0.5 K. The behaviour 
is reminiscent of earlier measurements on Bi-2212 7], 
also shown as an isotherm at 6 K, where the high-field 
"plateau" is preceded by a "kink" . Such similarity, mea- 
sured here for the first time in YBCO, is suggestive of a 
common origin. Since the present work is on a sample of 
the highest purity and in field-cooled measurements, this 
phenomena cannot simply be dismissed as material de- 
pendendent extrinsic behaviour, nor solely related to gap 
anisotropy [23|. The hysteresis seen in Bi-2212 measure- 
ments §\ is also naturally explained by the dominance of 
vortex scattering. The dramatic difference in magnitude 
between these two measurements is likely a consequence 
of the small electronic conductivity relative to the phonon 
contribution (which has not been subtracted) at 6 K in 



Bi-2212 and the comparatively huge electronic contribu- 
tion measured in this high-purity YBCO. 

An alternative explanation has recently been offered by 
Franz and Vafek [23| . In their fully quantum-mechanical 
theory, the Meissner state (at zero field) and the vortex 
state emerge as two distinct c?-wave states with different 
quasiparticle effective velocities. They both exhibit uni- 
versal conductivity, with different values of the universal 
limit (see Fig.|21). This appealingly accounts for the fact 
that the conductivity of the two different YBCO samples 
is the same not only at zero field but also at high fields 
(see Fig. It is not clear, however, why the finite tem- 
perature correction to this universal limit should be so 
dramatically different in the two states. 

In conclusion, when analyzed in terms of the quasiclas- 
sical theory of transport for a d-wave superconductor, 
the thermal conductivity of an extremely pure sample 
of YBCO reveals two features: 1) the universal limit as 
T ^ is confirmed, 2) the usual assumption that im- 
purity scattering can be treated as single isotropic phase 
shift in the unitary limit (or the Born limit) is incorrect. 
Transport theory as it stands must be revised, at least in 
the clean limit, perhaps by going to intermediate phase 
shifts and maybe in more profound ways. Moreover, in 
the presence of a magnetic field, we find that transport 
appears to be rapidly dominated by vortex scattering, 
which can therefore not be neglected as it usually is. 
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